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Abstract 

It is known that the separation of electrons into spinons and chargons, the spin-charge separation, 
plays a decisive role when describing strongly correlated one-dimensional (ID) Friedel oscillations. 
Here, we extend the investigation by considering a third electron fractionalization: the separation 
into spinons and orbitons. Specifically, we deal with two exact constraints of exchange-correlation 
(XC) density-functionals: (i) The constancy of the highest occupied Kohn-Sham eigenvalues upon 
fractional electron numbers, and (ii) their discontinuities at integers. By means of ID Hubbard 
chains, we show that spin-orbital separation can be decisive when dealing with derivative discon- 
tinuities of XC potentials, especially at strong correlations. 

PACS numbers: 71.10.Fd, 71.15.Mb 
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I. INTRODUCTION 



Spin and charge use to be treated as fundamental properties of ordinary electrons. How- 
ever, as predicted in the 1960s, when confined in one dimension interacting electrons display 
the unusual property of separating their spin and charge into two independent quasiparticles: 
spinous and chargons.- Both behave just like ordinary electrons, but: spinous have spin-1/2 
and no electrical charge, while chargons are spinless charged electrons. In the 1980s,-i^ an 
additional fractionalization was shown to occur: the spin-orbital separation, for which spin 
and orbital degrees of freedom are decoupled to form the orbitons - particles with no spin and 
charge, carrying solely the orbital information. The spin-charge and spin-orbital separations 
have recent evidences of experimental observation,-!^ teaching us that ordinary electrons can 
be considered bounded states of spinous, chargons and orbitons, three quasiparticles which 
can move independently and assume different velocities. 

The Kohn-Sham (KS) formalism of density-functional theory (DFT),-i^ by construction, 
retains the spin, charge and orbital degrees of freedom together, once it considers an aux- 
iliary system of noninteracting particles. It has been shown that the separation into char- 
gons and spinous are decisive when dealing with strongly correlated one- dimensional Friedel 
oscillations.- In this context, W. Yang et. al. have assumed strongly correlated systems as 
one of the modern challenges for DFT: "The challenge of strongly correlated systems is a 
very important frontier for DFT. To fully understand the importance of these systems, they 
must he looked at within the wider realm of electronic structure methods. Here, we intend 
to give a contribution into this challenge: By means of the one-dimensional Hubbard model 
(IDHM),- we investigate the effect of electronic correlation on the spin-charge and spin- 
orbital separations. Specifically, we consider two exact constraints of exchange-correlation 
(XC) density-functionals: (i) the constancy of the highest occupied Kohn-Sham eigenvalues 
upon fractional electron numbers, and (ii) their discontinuities at integers. These constraints 
are usually not satisfied even by modern approaches, and are the cause of dramatic errors 
when dealing with any generic situation involving transport of charges. For example, the 
absence of a proper derivative discontinuity leads local and semi-local functionals to fail in 
yielding neutral atoms in molecular dissociation processes: instead, they predict dissociated 
subsystems with unphysical fractional charges.— This is usually associated to a difficulty 
in producing localized orbitals, resulting in excessive charge transfer which is sometimes 
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referred to as delocalization error— of density functionals. 

In detail, we shall compare the performance of the IDHM local-density functional, its 
Perdew-Zunger self-interaction correction (PZSIC) and its spin-charge separation correction 
(SCSC). At strong correlations, we show that spin-orbital separation can be decisive when 
dealing with derivative discontinuities of XC potentials. 

II. THEORETICAL BACKGROUND 

(1) The one- dimensional Hubbard model: In one dimension, and in second-quantized nota- 
tion, the Hubbard model^ is defined as 

L L L 

H = -tJ2 (4^9+1,- + H.c.) +UJ2 4t^nc]iCji + J2 (1) 
i.c j j,o" 

where L is the number of sites, t is the amplitude for hopping between neighboring sites and 
Vj is a local external potential acting on site j. Occupation of each site is limited to two 
particles, necessarily of opposite spin. For this reason, in spite of using the XC nomenclature, 
we shall deal only with correlation, especially in situations for which its major counterpart 
is due to the local (on-site) interaction U. 

The basic Hohenberg-Kohn and KS theorems of DFT hold for the IDHM too, once the 
density n{r) is replaced by the on-site occupation n(j).^^ In terms of this variable, local- 
(spin)-density approximations for Hubbard chains and rings have been constructed.—"— 
Here, we chose the BALDA/FN^^ as the reference functional, considering only non-polarized 
systems with = N^. 

(2) Self -interaction correction: The BALDA/FN suffers from a self-interaction error. Perdew 
and Zunger have proposed the most known and used self- interaction correction (SIC),— 
written as: 

o-=t,4- k 

where nka{r) = /fco-|'?/'fco-(r)p, for the KS orbitals ipka and with occupation numbers 
fkcr- The approach of Perdew and Zunger has a long history, including statioii'^^'^ and 
time-dependent calculations.— Once the PZSIC functional of Eq. ([2]) depends explicitly 
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on the KS orbitals, via the orbital densities nka, we here shall employ the optimized 
effective potential (OEP) method^^ to obtain the orbital-independent PZSIC KS potential 
■^xc^CT^^NO) acting on each site j. 

(3) Fractional electron numbers: In a system with M = N + w{Q<w<l) electrons, the 
total ground-state energy is given by2^ 

E{N + w) = {l-w)E{N)+wE{N + 1). (3) 

Assuming that only the highest occupied (HO) KS orbital can be fractionally occupied, 
Janak^^ has proved that 

5E 

—— = €ho{N) = constant, (4) 
oM 

where Eka are the KS eigenvalues. The derivative discontinuities of the energy at each integer 
is then given by: 

E,{N) = IP{N) - EA{N) = BM+iiN + 1) - £^(iV), (5) 
which is also referred to as fundamental energy gap.— 

(4) The spin-charge separation correction: Considering spinous and chargons as independent 
entities, the occupied states will follow the semion distribution^^ of Fig. [I](b), where charge 
and spin were fractionalized to form the spin— 1/2 spinous (as a sum of two 1/4 spins) and 
the spinless chargons (as a sum of two — e/2 charges). As shown in Fig. [U^c), it has been 
proposed^, that the occupied states of a non-interacting KS system are built by retaining 
spin and charge together, at expense of the presence of holons (the chargon antiparticles) , 
whose densities are given by p^(j). The SCSC is then written as: 

-^^[p1(j)-^."c^"°^[p1(j) 

= VextU) + VH[n]i3) + t;L^'^M(j), (6) 



with 

N 



2 

a=t,i 



W 



(7) 
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(b) Interacting system (c) Kolin-Sliam system 




spins + cliarges spinons + chargons spins + cliarges togetlier 
togetlier + orbiton separated + spinless liolons + A 



i j i=> fractional occupation 

FIG. 1: (Color online) Schematic representation of occupied orbitals for 4 < < 5. Fractional 
occupations are indicated by the dashed squares. Panel (b): The orbiton is indicated by the orbital 
excitation from level k = 3 to k = 4. 

and 

cr=t,4. k 

where < w < 1, allowing fractional occupation. 
III. RESULTS 

Considering open Hubbard chains, in Figs. [2] (a)-(c) we plot euo versus fractional charge 
occupations for three values of U. We note that BALDA and PZSIC yield a linear depen- 
dence for enO) indicating a parabolic behavior for the total energy functional. The slope 
of each straight line tend to increase as U is increased. In addition, while the exact (on 
Lanczos precision) energy gaps become larger as U is increased, the values of Eg yielded by 
LDA and PZSIC remain almost the same for all values of U, scaling to ~ 50% of the exact 
values when U ^ l.Ot. The SCSC, which includes the spin-charge separation mechanism 
and is also free from one-electron self-interaction error, performs much better than both, 
yielding almost constant values for eno; even though the values of Eg remain also constant 
when increasing U. 

This difficulty in yielding accurate energy gaps is intrinsic of most available density func- 
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tionals. In this context, based on the schematic representation of Figs. [U^b) and[T]^c), we 
propose an alternative approach. Beyond spinons and chargons, in Fig. [^b) we indicate an 
excitation from level k = 3 to k = 4. This is the schematic representation of an orbiton: an 
orbital excitation which behaves like a spinless and uncharged electron. When dealing with 
the noninteracting KS electrons of Fig. [U^c), the energy gaps should then be increased by a 
constant A, which is equivalent to include the presence of an "anti-orbiton" , the extension 
of the SCSC idea of Eq. IQ. Thus, in a system with spinons + chargons + orbitons, we 
propose the XC potential to be given by 

approx+A^^^ (j) = <f (j) + A, (9) 



with 

As conceptually thought, noninteracting systems should have f] = 0, indicating the absence 



A = 4S r/ = 0,l,2,... (10) 



of spin-orbital separation. Strongly correlated systems, on the other hand, should have 
r]>l. 

In Figs. M (d)-(f) we plot the results for the BALDA+A, PZSIC+A and SCSC+A ap- 
proaches. When dealing with U < 3t, the best choice is 77 = 0. For stronger interactions, 
the values of Eg given hj rj = 1 are in great agreement with the exact values, multiplying 
the results of Figs. [2] (a)-(c) by a factor of ~ 2. Possibly, for U ^ lOt, it can be necessary 
to use 1] > 2. This is a clear delimitation between weakly and strongly correlated systems, 
whose connection between the increament of both, correlation and XC discontinuities, is in 
accordance with other well-known examples: In dissociation processes, to preserve neutral- 
ity of isolated atoms, the increment of atom-atom separation uses to be followed by a joint 
increment of electronic correlation and XC discontinuity.- 

In Figs. [3] (a)-(c) we plot the XC potential profiles for U = 5.0t and two slightly different 
electron numbers. When dealing with strong interactions, it has been shown that the XC 
potential should oscillate with Akp-— Here, only SCSC yields the correct frequency, as seen 
by the four negative peaks in Fig. [3]^c). The orbiton correction, as given by Eq. (|9]) with 
?7 = 1, induces a correct constant discontinuity in the XC potential when crossing = 4, 
as seen by the jump from A^ = 3.999 to A^ = 4.001. The same is not observed for the XC 
potentials with = 0, for which the curves are almost identical. 
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Density profiles: The density profiles of ID Hubbard cliains, wliicli are invariant under 
tlie correction of Eq. (|9]), are known to suffer a 2/ci? — )■ 4kp crossover as U is increased.- 
Here, maintaining U fixed, we investigate tlie density profiles upon variation of electron 
numbers. In Figs. |3] (d)-(f) we show the density profiles for 2 < < 4, with an increase of 
0.2 electrons for each curve, from bottom to top. The dashed lines indicate exact values (on 
Lanczos precision) for N = 2 and = 4. The 4k p Wigner crystal oscillations (identified by 
two or four peaks) are consequence of the strong interaction U = 10. Ot. Only SCSC is able 
to reproduce the correct density profiles. Note that a similar change in frequency (from two 
to four peaks) could be obtained by fixing = 4 and varying the interaction from f/ = to 
t/ — )■ oo.- Similar conclusions remain valid when including a confining harmonic potential 
(not showed here). 

The emergence of Akp Wigner crystal oscillations, via Eq. (Q, seems to be connected 
with the necessity of using Eq. (Q, with rj > 1. Both situations have a common point: The 
strongly interacting regime. W. Yang el. al. have recently pointed out^ that, beyond strong 
correlations, it is a basic challenge to understand whether the KS orbitals and eigenvalues 
have any further significance. Therefore, the spinon-chargon-orbiton separation can be a 
first step into this direction. 

In summary, specifically for the BALDA functional, we can conclude that: (I) spin-orbital 
separation, when included by considering the orbiton excitations, yields an increment of en- 
ergy gaps; (II) spin-charge separation, when included by means of the SCSC XC potential, 
yields ikp density oscillations and almost constant highest occupied KS eigenvalues. Consid- 
ering (I) and (II), we can argue that when dealing with strong correlations, electrons should 
not be treated as unique particles. Instead, the separation into spinous, chargons and or- 
bitons must be considered. The way to do so in a noninteracting KS calculation, which by 
construction retains spin, charge and orbital degrees of freedom together, is a combination 
of (I) and (II). 

IV. CONCLUSIONS 

Constancy of eho- The BALDA and BALDA+PZSIC yielded typical results attributed to 
delocalization errors of density functionals: Incorrect linear behavior of ^ho upon fractional 
charge occupation, with an increase in the slope when interaction is increased. On the 



7 



CO 



0.4 
0.3 
0.2 

O 

^ 0.1 
0.0 
-0.1 



L=13 U=1.0t 



Exact 

BALDA 
X PZSIC 
A SCSC 



3.0 



3.5 



4.0 

N 



4.5 



5.0 



+ 

C 

to' 



0.4 
0.3 

0.2 

oO.i 

0.0 
-0.1 



(d) L=13 U=1.0t 



Exact r 

BALDA + A I 
X PZSIC + A I 
A SCSC + A I 



3.0 



3.5 



4.0 

N 



4.5 



5.0 



0.5 
0.4 
0.3 
0.2 
X 0.1 
0.0 
-0.1 
-0.2 



CO 



(b) 



L=13 U=5.0t 



3.0 



3.5 



4.0 

N 



4.5 



5.0 






0.5- 




0.4- 




0.3- 


< 


0.2- 


+ 


0.1- 


o 




CO 


0.0- 




-0.1- 




-0.2- 




0.6-| 




0.5- 




0.4- 




0.3- 


< 


0.2- 


+ 


0.1- 


o 




X 

CO 


0.0- 




-0.1- 




3.0 



-0.2 




FIG. 2: (Color online) Highest occupied KS eigenvalues eho obtained by means of Eq. ([9]). Panels 
(a)-(c): T] = 0; Panels (d)-(f): r] = 1. 
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FIG. 3: (Color online) (a)-(c): XC potential profiles for two slightly different electron numbers. 
(d)-(f): Density profiles for 2 < < 4. The dashed lines indicate exact values for N = 2 and 
= 4, from bottom to top. 
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other hand, the BALDA+SCSC, which has been especially conceived to deal with strong 
interactions, yields almost constant values for eno- The separation into spinons + chargons 
seems to be more efficient in correcting the BALDA than the Perdew-Zunger prescription. 

Derivative discontinuity: Even though all functionals tested here yield accurate energy 
gaps at weak interactions, the accuracy is lost when interaction is increased. By means 
of Eq. (Q, we proposed to also consider orbiton excitations. The results indicate a clear 
delimitation between weakly and strongly correlated systems: At strong correlations and 
with A 7^ 0, the inclusion of spin-charge and spin-orbital separations yields accurate energy 
gaps with correct 4kp density proffies. It remains to be seen whether the same trends are 
observed in ab initio calculations under non-local 1/r Coulomb interactions. Understanding 
possible connections between using A 7^ and other types of discontinuity corrections is a 
topic for further investigations. 
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